EFI 08-27 



Constraining Light Bosons with Radiative T(IS') Decays 

David McKeeifl 

Enrico Fermi Institute and Department of Physics, 
University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637 
(Dated: March 30, 2009) 

Light bosons can be found in large classes of theories beyond the standard model. These light 
bosons may not be ruled out by current experiments and, indeed, may even provide an explanation 
for some anomalous experimental results. The radiative decays of quarkonium (cc, bb) states offer 
a promising opportunity to investigate such light bosons. Specifically, we investigate the reach that 
current CLEO data can have on models with light scalar and pseudoscalar bosons. 
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^ ; I. INTRODUCTION 

I Models that contain light CP-even and CP-odd bosons occur in many classes of theories beyond the standard model. 
' They are a generic feature of models with an expanded Higgs sector. The extended freedom in the scalar mass matrices 
' allows a wide range of parameter space in which some bosons may be light while reducing their couplings to the Z 
. boson so that LEP would not have observed them. A particularly attractive example is the next-to-minimal standard 
Qh' supersymmetric model (NMSSM) which solves the /i-problem of the MSSM. In the NMSSM, light pseudoscalars 
^ , can arise naturally [p. Minimal model- independent light dark matter scenarios also can involve light scalars or 
^ pseudoscalars as mediators to serve as an efficient annihilation channel. Such light bosons should mix with CP-even 
f~| or -odd Higgs bosons to couple to standard model fermions. Light scalars and pseudoscalars can also arise generically 
in models that contain light sgoldstinos 2] . There are also extra-dimensional gaugephobic models in which the Higgs 
is decoupled from gauge bosons, allowing it to be hght [1, |3|- 
^ . There are hints that light bosons may already have been observed experimentally. First, if a model contains a 
' pseudoscalar with mass, nia < ^rrii,, the LEP limit on the mass of the standard model Higgs boson {rrih > 114.4 GeV) 
OO 1 can be evaded. This could also provide an explanation for the slight excess of events seen at rrih — 100 GeV Q 
C ■ ' since the Higgs' branching ratio to b quarks is reduced !6]. In addition, these models are also able to accommodate 
' a light dark matter particle that could explain the excess of 511 keV photons from the Galactic center seen by the 
INTEGRAL satellite 7, ^ although recent measurements may cast the explanation in terms of dark matter in some 
doubt [9]. The DAM A/LIBRA Collaboration's observations JX)] can also be reconciled with null results from CDMS 
if dark matter is relatively light [Tl|j E^]- These scenarios are discussed further in Ref. [l3|. In addition, a model 
' containing a pseudoscalar boson with a mass of 214 MeV 14| can explain the three E+ pjX^ ^~ decays of the same 
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dimuon mass seen by the HyperCP Collaboration [15| . The probability of such a cluster of events in the standard 
model is regarded as below 1% [l^. Explanations of this signal as due to a light sgoldstino have also been studied [l3|. 
^ , However, recent results from CLEO rule out the explanation of the HyperCP results as being due to a pseudoscalar 
' of mass 214 MeV [l^. Furthermore, there are hints that light bosons could be the cause of a 3(7 deviation from the 
■ ■ ' SM prediction for the process 7r° e+e" at KTeV [li,[23,l2H. 

In this paper we investigate the implications of light scalars and pseudoscalars for the decays T(IS') 
7{/i+/i^, TT+TT^, if +iir^} and estimate the backgrounds from the radiative return processes e^e~ — *■ 77* 
7{/i"'"/i~, TT+TT", if'^ii'"}. The radiative return to a muon pair has recently been considered in Ref. ^25|. 

The case of light pseudoscalars is presented in Sec. II. The couplings of the pseudoscalar to standard model fermions 
and the expression for the signal they generate are given in Sec. HA. The background to this process is shown in Sec. 
IIB. In Section IIC, the experimental reach with respect to the parameters is discussed. This is repeated for the case 
of light scalars in Sees. IIIA, IIIB, and IIIC. In Sec. IV we conclude. 
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FIG. 1: Plots of Fa and Fs [Eq. ©] as functions of x = 2E-,/mr. 



II. PSEUDOSCALARS 
A. Signal (including scalars) 

We consider a model with a pseudoscalar A and a scalar S coupled to up- and down-type quarks and leptons. We 
assume the couplings to be flavor diagonal. We parameterize the interaction Lagrangian as 

Ca = ~ (^-^ [gdrndl^l + gdmad^^d -\- gumuUj^u) , (1) 

£5 = - (XdmeU + Xdrriddd + A„m„uM) . (2) 

with copies for each family. For simplicity, we have used equal coupling to leptons and down-type quarks, motivated 
by grand unified theories. The results should be straightforwardly generalized to different couplings. Here, v — 
(V2Gf)~^/^ ^ 246 GeV. FCNCs are generated by one loop two-quark operators and by four-quark operators. Limits 
from flavor changing decays such as K -nfi^ and B Kfj,~^fi~ will tend to provide stricter constraints than 
those from T decays. However, there are examples of models where parameters can be chosen so contributions to 
these processes can cancel. In particular, in the NMSSM, there exist ranges of squark mass splittings and chargino 
masses where such cancellations occur. For details, see Ref. For our purposes we will neglect flavor changing 

decays to concentrate as an independent probe on the limits that can be derived from T decays. The rate of the 
radiative decay of the T to a pseudoscalar is [1^ [l^l 



r(T ^ jA) _ gjmlGF f 
r(T^e+e-) V2™ 



Ca{.x) (3) 



and similarly for a scalar. 



r(T^e+e-) V27ra V m 



ix) (4) 



Here, Cp = [N^ ~ l)/2N = 4/3 for S'C/(3) and x 1 - m\ g/m\ = 2E^/mr. Ca{x) and Cs{x) describe QCD and 
relativistic corrections ,^24] and are given by 



Ca,s{x) = 



1-['^]Fa,s{-) 



(5) 



Fig. [T] shows Fa{x) and Fs{x). The differential rate is distributed with respect to the photon's center-of-mass angle 
with the beam as 1 -I- cos^ 6^ . 

For ruA < 2mr, we take B{A j.i'^ii^) ~ 1. The rate to two photons is loop suppressed and CP conservation 
forbids the s-wave decay to pairs of pseudoscalar mesons. Decays with three particles in the final state are taken to 
be negligible. 
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B. Background 

The dominant background to the process T(IS') — > ^fi'^ comes from the radiative return e+e^ — > 77* 

7/i+/x^. The differential cross section is ^2] 

da f '^'^ \ d cos 9^dxdip 
dml^ \ATrsJ {ml^){p^- ■ k){p^+ ■ k) 

+ ^ [{Pe- ■ (P^- +P^.+ ) f + {Pe+ ■ {P^.- 

with s the center-of-mass energy squared, k the four momentum of the photon, ip the azimuthal angle between the 
photon and muon, and x = 2E^- 1 \f^- This is plotted in Fig. [2] as a function of dimuon mass for several different 
limits on the photon's angle. 
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C. Probing the coupling to fermions 

One can set a lower limit on the value of B(T — > A'^')B{^A that is observable, or equivalently on g^, by 

equating the number of events due to the signal to the number due to statistical fluctuations of the background in 
a certain range of dimuon mass. We focus on the CLEO experiment, taking an integrated luminosity at the T(IS') 
resonance of 1.06 fb"^ HI], 21 X 10^ T(IS') produced, a resolution of 5 MeV on the dimuon mass, and a nominal 
efficiency of 2%. For different values of the pseudoscalar mass we plot the lower bound on the observable value of 
B(T 7a)S(A in Fig. [31 Indeed, CLEO reports limits that are roughly consistent with these estimates in 

Ref. [Tsl, in the range of (several) 10~^ over the mass range in Fig.[3l 

If we take B{A fj,~) ~ 1, this can simply be translated into a limit on the lowest possible value of gd observable 
by using Eq. ^ to express B{T jA) in terms of g^,- We use as{Mr) — 0.2 and note that, for ma,ms = 0, the 
QCD corrections given in Eq. (O are around 40% and 65% of the tree-level rates in the pseudoscalar and scalar cases, 
respectively. This indicates that there is a substantial amount of theoretical uncertainty and quantities obtained using 
Eqs. ^ and ^ could vary quite a bit. Figure [J] shows the limits on gd with the tree level QCD corrections taken 
into account. 

Figure [3] shows that there is a slight dependence on the cut on the photon's angle with the beamline and that an 
optimal value occurs somewhere in the range of 0.75 < cosOj < 0.85. This can be understood as being due to the 
angular distributions of the signal and background. The signal has a distribution proportional to 1 + cos^ 6^ while 
the background is roughly distributed as (1 — cos^ ^7)"^- If we denote the cut on cos 6'-y by c, the lower bound on gd 
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FIG. 3: Lower bound on observable values of S(T — > A^)B{A —> ^'^fi~) for a range of pseudoscalar masses with several values 
of cut on photon angle. It is seen that taking —0.85 < cosd^ < 0.85 generates a slightly stronger limit than either of the other 
two cuts. 
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FIG. 4: Lower bound on observable values of gd for a range of pseudoscalar masses with the photon cut near its optimum: 
-0.85 < cosO^ < 0.85. 



depends on c as 

This function has a maximum at c ~ 0.88. The background's deviation from the (1 — cos^ ^7)^^ distribution likely 
accounts for the minimum occurring at a shghtly lower value of c. 



III. SCALARS 



A. Signal (additional features) 

We now consider a model with a scalar S coupled to fermions as in Eq. [51 We define the ratio of the scalar 's up- 
and down- type couplings as r — Xu/Xd- The rate of the radiative decay of the T to a scalar is given by Eq. |4l As 
mentioned above, the differential rate is distributed as 1 + cos^ 9.y. 

In the case of the light scalar, determining its branching ratios is slightly more complicated than for a pseudoscalar. 
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The rate to muons is 

r(5-MV-) = ^^/3^ (8) 

with = (1 — 4m^/m|)^/^. We estimate the couphng to two pseudoscalar mesons following Refs. [1^ [l^l and 
references therein. Heavy quark (c, 5, t) loops generate a coupling of the scalar to two gluons 



S 



C.ff = j^\^-j{2Xu + Xd)G^^,G'^^'' . (9) 
The conformal anomaly then relates this term to the trace of the energy momentum tensor 

Q'; = ~^-^g;^g^'^''+ J2 ^^^^^^ ' (10) 

where 5 = 11 — 2np/i is the coefficient of the QCD (3 function, (3{gs) — —bg^/lGiT^. Putting this all together one 
obtains the interaction Lagrangian 

Ant = f ] ^J^h^ I ^ ^^^^^^ _ e;; I . (11) 



i i—u.d,s 



One can then use a chiral Lagrangian and make the replacement 

J2 mS^^^^^TT{^Hd)+h.c. (12) 

i—u^d,s 

with E a matrix containing the pseudoscalar octet, S = exp[{2i/ f.^)\°'TT°'], mo related to the quark condensate as 
(qiqi) ~ f^mo, and M a matrix containing the quark masses, M — diag(m„, m^, mj). Similarly we replace 

J2 A,m,Vi^^, ^^^Tr(EtM')+h.c. (13) 

i—u,d,s 

with il/' = diag(Aum„, Xdm^, Xdirig). We write the trace of the energy momentum tensor as 

= -2d''TT+d^n- + AmlTT+TT- + ... . (14) 



We then find 



(tt+tt I ^ mitlJi^i\0) = nio{mu + rud) 

i—u,d,s 

^rnl± , (15) 
(tt+tt"! X! AimiVJiV'ilO) = Admo(rm„ + md) 

i—u,d.s 

^^[{l + r)ml^ + {l-r)^ml] (16) 

where Ato^ = m^^o — m^^. (Electromagnetic shifts of the pseudoscalar masses have been ignored for simplicity.) 
Using Eq. (|14p one can evaluate the matrix element involving the trace of the energy momentum tensor as 

{7r+7T-\ef^\0)~q^ + 2ml±^ml + 2mli . (17) 
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FIG. 5: B{S — > yC^ as a function of ms for several values of r = A^/Ad, obtained from Eqs. (|18|l and (I20p . 



This is used to find 



r(5 ^ TTvr) mi^ V V 4 

+ ^(1_,) 



4 mi; } 

with /^TT = (1 — 4m^/m|)^/^ and both charged and neutral pions included. If we set r = 1 wc find 



T{S ^ tttt) mi; { V 2 / m;| J (3^ 



(19) 



which agrees with Eq. (2.50) of [26]. However, above kaon threshold the chiral Lagrangian picture may no longer hold. 
For ms > 2mK we follow Ref. [26] and use a perturbative spectator model to estimate the scalar's partial widths. 
Denoting T{s /i+/i^) as F^^, etc., we take 

'■ ^uuMd ■ ^ss '■ Tgg : ^NN = ^dm^^f^f^ : 3A^j_^m^ ,^/?^ (20) 

: --{^y^ + A^) - (A^ + X^Pl ~ A^/?^} : SA^ .m^ ,/3^ 

where N denotes a nucleon. We equate the amplitude for pio n production to that of nucleon production, ignoring 
differences between diquark and quark pair creation. As in [26|, we determine m„, m^, ms, and as by matching the 
chiral model onto the perturbative spectator model in their overlapping regions of applicability. We find that to,; = 
md = 50 MeV, ms = 450 MeV and as — 0.47 allow the perturbative spectator model to reproduce the chiral model 
fairly well. Current quark masses do not allow for as good agreement between the two models. We ignore the u and d 
quark content in the kaons and the 77 and approximate F^+^- +F^o^o ~ F„g + F^j, Tj^+j^- -j-Tj^oj^o +F^^ ~ ^ss- We 
use flavor SU{3) to estimate T^+j^- ~ 2F^o^o and Tj^+j^- ~ Fj^oj^o ~ (9/8)F^^ above rj threshold. Note that we have 
used the wavefunction ry — {uu + dd + 2ss) /\/6. Using the more phenomenologically motivated 7] ~ {uu + dd + ss)/ V3 
will decrease the importance of the 7777 channel relative to kaons. For simplicity, we ignore tt'^-t] production, which 
vanishes in the limit A„mu — Xdm^. Higher multiplicity final states should not enter into the picture until p pair 
production threshold is crossed. We assume that these final states are phase-space suppressed. In this approximation, 
we find that the scalar's branching ratio to baryons is always less than 3%. We include baryon production in our 
determination of other branching ratios but do not plot it for clarity. Further studies with these assumptions relaxed 
are warranted. We plot the estimates of the branching ratios of the scalar to muons, pions and kaons as a function of 
ms for several values of r in Figs. OHl 
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FIG. 9: Scalar's branching ratios as a function of ms for r = \u/\d = 1, obtained from Eqs. f|18|l and f|20|l . Baryon production 
is included in this calculation but is found to be less than three percent for all values of ms and is not plotted for clarity. 
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FIG. 10: The cross section for e^e~ 77* "f-K^n" with —0.85 < cosS-y < 0.85. 



B. Background 

We take as the dominant background for T(IS') — > 76" — > 77r+7r^ the radiative return process e^e^ — > 77* ^ir^-n^ 
and similarly for the background to T(IS') — > 76" — > jK^K" . We use vector meson dominance to write the pion and 
kaon electromagnetic form factors as 

with Npyr = 1, N^yr = N^^ = 0, NpK = 1/2, N^^x = 1/6, and N^^ — 1/3. The differential cross sections, da/dm^j^ 
and da/diriKK, are plotted in Figs. [TUl and [TTl 



C. Probing the coupling to fermions 

We can perform a similar analysis of the radiative return background to pion and kaon pairs to set a lower limit 
on the observable values of B{T jS)B{S tt+tt") and ^(T ^ jS)B{S K+R-). These are plotted in Figs. [H 
and [T31 We can then use Eq. ^ with the branching ratios of S obtained from Eqs. and (|^D|) to translate the 
branching ratio limits to limits on observable values of Ad- We use the same naive efficiency and resolution as in the 
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FIG. 11: The cross section for e+e" -> 77* -fR-^K' with -0.85 < cos (9^ < 0.85. 
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FIG. 12: Lower bound on observable values of B{T Sy)B{S n^n~) for a range of scalar masses with —0.85 < cos^^ < 0.85. 
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FIG. 13: Lower bound on observable values of B(T S"f)B{S K^K~) for a range of scalar masses with —0.85 < cos 9^ < 
0.85. 
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FIG. 14: Lower bound on observable values of obtained from the final states indicated with —0.85 < cos 6^ < 0.85 for r • 
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FIG. 15: Lower bound on observable values of Ad for r — 0.1. 



pseudoscalar analysis. The limits that come from 7/1+/^", 777 "•"tt", and jK~^K~ final states will depend on the value 
of r = Xu/^d- These limits, for several values of r, are plotted in Figs. [TUfTCl We observe that the hmit due to the 
pion final state becomes more stringent as r is increased. For r > 1 the dimuon final state does not offer very useful 
bounds and the hadronic final states become more important. 

For large m^^ or tukk, the background is suppressed since the electromagnetic form factors are far from resonance. 
Thus, we see that we can probe down to smaller couplings than in the dimuon case where the increased resolution 
becomes a limiting factor. At large enough mTTTr or ttikki we do not expect any background. Requiring that we see 
three signal events limits A^ > 0.2. This can be thought of as an absolute lower bound on observable values of the 
coupling. 

Assuming that the dominant background is due to the radiative return process, the minimum coupling that can be 
probed in both the scalar and pseudoscalar cases scales as 

(l laJCdtV 

where N is the number of quarkonia produced, cr is the cross section of the radiative return background, J £ dt is the 
integrated luminosity, and e is the efficiency. At a given q^, a ^ s^^ = Aty'^, where My is the mass of the resonance. 
Then the minimum coupling scales as {eMyN"^ / J C dt)~^^^. 
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FIG. 16: Lower bound on observable values of for r = 10. 



IV. CONCLUSIONS 



Radiative decays of the T(IS') can provide useful information about light scalars and pseudoscalars. There is 
solid motivation for considering light spinless particles. A model-independent extension of the scalar sector can 
generically include light scalars or pseudoscalars that would have escaped detection due to the increased freedom 
in the couplings of the scalar sector. Such models can accommodate light dark matter that may be present. Such 
scalars also occur in the NMSSM, models with light sgoldstinos, and "gaugephobic" Higgs scenarios. In addition, 
light spinless particles may help to shed light on several anomalous experimental observations. Quarkonium offers 
an excellent environment to probe the tree level coupling of light scalars and pseudoscalars to fcrmions and provides 
results that are complementary to those that involve flavor changing couplings. These results can help to tease out 
the nature of a scalar sector beyond the standard model. 

We have seen that current CLEO data on the T(IS') resonance can limit couplings down to gd,^d > 0.25 in the 
Lagrangians of Eqs. ([T]) and In particular, they are able to rule out [l3| the interpretation of the HyperCP signal 
as due to a light pseudoscalar Higgs which requires ~ 0.4 (l3|- Looking at Fig. [TH current CLEO data can also 
directly rule out a standard model Higgs boson (A^ = 1, r = 1) up to masses of about 1 GeV in the dimuon channel. 
This is subject to some uncertainty in the branching ratios, although chiral perturbation theory should not be suspect 
in this area. 

Other models can also be constrained. In models with light sgoldstinos the coupling in the interaction Lagrangian 
of Eq. [T]is given hy g/v — Aj/F and in Eq. [Uby X/v = Aj /F where y/F is the scale of SUSY breaking and Aj is 
a soft trilinear coupling. Limiting gd, Xd > 0.25 constrains F/Af > 1 TeV. If one sets Af ~ 100 GeV then one can 

limit \/F> 315 GeV. _ 

In gaugephobic Models, one typically finds < 0.1 [4]. Current data are seen above to approach this limit. Further 
study would be welcome. 

Stated more generally, one can imagine the interactions of Eqs. (jlj and ^ suppressed by the scales — v/gd, 
As — v/Xd which are probed up to Ka.s — 1 TeV. The scale that can be reached at another resonance should go as 
^A,s ~ {^MyN^/ J C dtY/^. The theoretical uncertainty involved in obtaining these results is fairly large, however, 
and findings could differ somewhat from those presented here. In particular, a more rigorous incorporation of hadronic 
decays could help shore up conclusions that can be drawn about light spinless particles coupled to fermions. 
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